Plasmonics in graphene at infra-red frequencies 
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We point out that plasmons in doped graphene simultaneously enable low-losses and significant 
wave localization for frequencies below that of the optical phonon branch huioph ~ 0.2 eV. Large 
plasmon losses occur in the interband regime (via excitation of electron- hole pairs) , which can be 
pushed towards higher frequencies for higher doping values. For sufficiently large dopings, there is a 
bandwidth of frequencies from cooph up to the interband threshold, where a plasmon decay channel 
via emission of an optical phonon together with an electron-hole pair is nonegligible. The calculation 
of losses is performed within the framework of a random-phase approximation and number conserv- 
ing relaxation-time approximation. The measured DC relaxation-time serves as an input parameter 
characterizing collisions with impurities, whereas the contribution from optical phonons is estimated 
from the influence of the electron-phonon coupling on the optical conductivity. Optical properties of 
plasmons in graphene are in many relevant aspects similar to optical properties of surface plasmons 
propagating on dielectric-metal interface, which have been drawing a lot of interest lately because 
of their importance for nanophotonics. Therefore, the fact that plasmons in graphene could have 
low losses for certain frequencies makes them potentially interesting for nanophotonic applications. 

PACS numbers: 73.20.Mf,73.25.+i 



I. INTRODUCTION 



In recent years, an enormous interest has been sur- 
rounding the field of plasmonics, because of the variety 
of tremendously exciting and novel phenomena it could 
enable. On one hand, plasmonics seems to be the only vi- 
able path towards realization of nanophotonics: control 
of light at scales substantially smaller than the wave- 
length 0, H, H, 3- On the other hand, plasmonics is 
a crucial ingredient for implementation of most meta- 
materials, and thereby all the exciting phenomena that 
they support [E IE 0, 111 > including negative refraction, 
supcrlcnsing, and cloaking. However, there is one large 
and so far insurmountable obstacle towards achieving this 
great vision: plasmonic materials (most notably metals) 
have enormous losses in the frequency regimes of interest. 
This greatly motivates us to explore plasmons and their 
losses in a newly available material with unique proper- 
ties: graphene 0, EL EL 01 HI EL EI ■ 

Graphene is a single two-dimensional (2D) plane of car- 
bon atoms arranged in a honeycomb lattice, which has 
only recently been demonstrated in high qua lity samples 
and with superior mobilities 0, El, EL EL EL EL El- 
This material is a zero-gap semiconductor, which can be 
doped to high values of electron or hole concentrations by 
applying voltage externally @, much like in field effect 
transistors (FET). While this kind of control over electri- 
cal properties of materials is at the heart of modern clec- 
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tronics, it was also demonstrated that the same proce- 
dure (electric gating) [IE E} leads to a dramatic change 
in optical properties of graphene because of its impact on 
the strong interband transitions. Collective excitations 
(plasmons) in graphene ho ld p otential for technolo gica l 
applications as well [H, EL EL EL El, EL EL HE ELII3 ; 
for example, coherent terahertz sources based on plas- 
mon amplification were suggested and discussed in Rcfs. 
[l9L E| • Graphene was predicted to support a transverse 
electric (TE) mode [23, which is not present in usual 
2D systems with parabolic electron dispersion. Thcrmo- 
plas ma p olaritons in graphene have been discussed in 
Ref. [la ], pointing out at new opportunities in the field 
of plasmonics. 

Here we investigate plasmons in doped graphene and 
demonstrate that they simultaneously enable low-losses 
and significant wave localization for frequencies of the 
light smaller than the optical phonon frequency hwoph ~ 
0.2 eV 28[. Interband losses via emission of electron- 
hole pairs (1 st order process) can be blocked by suffi- 
ciently increasing the doping level, which pushes the in- 
terband threshold frequency uii nter towards higher values 
(already experimentally achieved doping levels can push 
it even up to near infrared frequencies). The plasmon 
decay channel via emission of an optical phonon together 
with an electron-hole pair (2 nd order process) is inac- 
tive for to < cooph (due to energy conservation), how- 
ever, for frequencies larger than iooph this decay chan- 
nel is nonnegligible. This is particularly important for 
large enough doping values when the interband threshold 
Winter is above Lu 0p h- in the interval uj 0p h < u < Winter 
the 1 st order process is suppressed, but the phonon decay 
channel is open. In this article, the calculation of losses is 
performed within the framework of a random-phase ap- 
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FIG. 1: (color online) (a) Schematic description of a surface plasmon (SP) on metal-dielectric interface, (b) SP dispersion 
curve (solid blue line) for Ag-Si interfaces; dotted blue is the light line in Si; dashed red line denotes the SP resonance, (c) 
Wave localization and propagation length for SPs at Ag-Si interface (experimental Ag losses are taken into account). 



proximation (RPA) and number conserving relaxation- 
time approximation p9l |; the measured DC relaxation- 
time from Ref. Q serves as an input parameter char- 
acterizing collisions with impurities, whereas the optical 
phonon relaxation times are estimated from the influence 
of the electron-phonon coupling [3(| on the optical con- 
ductivity (3lj |. 

In Sec. UH we provide a brief review of conventional 
surface plasmons and their relevance for nanophotonics. 
In Sec. Mil we discuss the trade off between plasmon 
losses and wave localization in doped graphene, as well 
as the optical properties of these plasmons. We conclude 
and provide an outlook in Sec. IIVI 

II. SURFACE PLASMONS 



lengths, which arc conveniently quantified (in terms of 
the SP wavelength) with the ratio ^R.q sp /Qq sp ; this quan- 
tity is a measure of how many SP wavelengths can an 
SP propagate before it loses most of its energy. The 
wave localization (or wave "shrinkage") is quantified as 
Aair/^sp, where \ a i r = ^c/uj (the wavelength in air). 
These quantities are plotted in Fig. [TJc) for the case of 
Ag-Si interface, by using experimental data (see Q and 
references therein) to model silver (metal with the low- 
est losses for the frequencies of interest). Near the SP 
resonance, wave localization reaches its peak; however, 
losses are very high there resulting in a small propaga- 
tion length I ?s 0.1A sp ~ 5nm. At higher wavelengths one 
can achieve low losses but at the expense of poor wave 
localization. 



Surface plasmons (SPs) are electromagnetic (EM) 
waves that propagate along the boundary surface of a 
metal and a dielectric [see Fig. ffl^a)]; these are trans- 
verse magnetic (TM) modes accompanied by collective 
oscillations of surface charges, which decay exponentially 
in the transverse directions (see, e.g., Refs. P, and 
Refs. therein). Their dispersion curve is given by: 

c y € r + e(w) 

[see Fig. D^b)]; note that close to the SP resonance 
(uj = ujsp), the SP wave vector [solid blue line in Fig. 
QJb)] is much larger than the wave vector of the same 
frequency excitation in the bulk dielectric [dotted blue 
line in Fig. |Tfb)] . As a result, a localized SP wave 
packet can be much smaller than a same frequency wave 
packet in a dielectric. Moreover, this shrinkage is ac- 
companied by a large transverse localization of the plas- 
monic modes. These features are considered very promis- 
ing for enabling nano-photonics [H, 0, [!, 3 , as well as high 
field localization and enhancement. A necessary condi- 
tion for the existence of SPs is e(uj) < — e r (i.e., t(uj) is 
negative), which is why metals are usually used. How- 
ever, SPs in metals are known to have small propagation 



III. PLASMONS AND THEIR LOSSES IN 
DOPED GRAPHENE 

Graphene behaves as an essentially 2D electronic sys- 
tem. In the absence of doping, conduction and valence 
bands meet at a point (called Dirac point) which is also 
the position of the Fermi energy. The band structure, 
calculated in the tight binding approximation is shown 
in Fig. 2(b) (see Ref. [25| and references therein); for 
low energies the dispersion around the Dirac point can 
be expressed as E n ^ = nvpU\\<i\, where the Fermi veloc- 
ity is vp = 10 6 m/s, n = 1 for conduction, and n = —1 for 
the valence band. Recent experiments [32j have shown 
that this linear dispersion relation is still valid even up to 
the energies (frequencies) of visible light, which includes 
the regime we are interested in. 

Here we consider TM modes in geometry depicted in 
Fig. 2(a), where graphene is surrounded with dielectrics 
of constants e rl and e r i- Throughout the paper, for defi- 
niteness we use e rl = 4 corresponding to Si02 substrate, 
and e r 2 = 1 for air on top of graphene, which corresponds 
to a typical experimental setup. TM modes are found by 
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FIG. 2: (color online) (a) Schematic of the graphene system and TM plasmon modes. Note that the profile of the fields looks 
the same as the fields of an SP [Fig. HJa)]. (b) Electronic band structure of graphene; to indicate the vertical scale we show 
the Fermi energy level for the case Ef = 1 eV. (c) Sketch of the intraband (green arrows) and interband (red arrows) single 
particle excitations that can lead to large losses; these losses can be avoided by implementing a sufficiently high doping, (d) 
Plasmon RPA and semiclassical dispersion curves. Black solid (RPA) and black dot-dashed (semiclassical) lines correspond to 
e r j — e r 2 = 1; Blue dashed (RPA) and blue dotted (semiclassical) lines correspond to e r i — 4 and e r 2 = 1. The green (lower) 
and rose (upper) shaded areas represent regimes of intraband and interband excitations, respectively. 



assuming that the electric field has the form 

E z = Ae iqz ~ QlX , E y =0,E X = Be lqz ~ QlX , for x > 0, 

E, = Ce lqz+Q2X , E y =0,E X = De lqz+Q2X , for x < 0. 

(2) 

After inserting this ansatz into Maxwells equations and 
matching the boundary conditions [which include the 
conductance of the 2D graphene layer, <r(uj,q)], we ob- 
tain the dispersion relation for TM modes: 



<72 



toe 



(3) 



By explicitly writing the dependence of the conductiv- 
ity on the wave vector q we allow for the possibility of 
nonlocal effects, where the mean free path of electrons 
can be smaller than g _1 [33|. Throughout this work we 
consider the nonrctarded regime (q 3> w/c), so equation 
© simplifies to 



q « Qi « Q2 ~ £o 



Erl + (r2 ZlUJ 

2 a(ui,q)' 



(4) 



Note that a small wavelength (large q) leads to a high 
transversal localization of the modes, which are also ac- 
companied by a collective surface charge oscillation, sim- 
ilar to SPs in metals; however, it should be understood 
that, in contrast to SPs, here we deal with 2D collective 
excitations, i.e. plasmons. We note that even though 
field profiles of plasmons in graphene and SPs in metals 
look the same, these two systems are qualitatively dif- 
ferent since electrons in graphene are essentially frozen 
in the transverse dimension (34[. This fact and the dif- 
ferences in electronic dispersions (linear Dirac cones vs. 
usual parabolic) lead to qualitatively different dispersions 
of TM modes in these two systems [see Fig. [IJb) and 
Fig. [DJd)]. To find dispersion of plasmons in graphene 



we need the conductivity of graphene a(u>,q), which we 
now proceed to analyze by employing the semiclassical 
model [33| (in subsection IIII A[) . RPA and number con- 
serving relaxation-time approximation [2t| (in subsection 
IIIIBI) . and by estimating the relaxation-time due to the 
influence of electron-phonon coupling [3fjj on the optical 
conductivity [3l[ (in subsection IIII C[) . 



A. Semiclassical model 

For the sake of the clarity of the presentation, we first 
note that by employing a s imp le semi-classical model for 
the conductivity (see Ref. [33[), one obtains a Drude-likc 
expression [35| : 



e 2 E F 



ttH 2 UJ + IT 



(5) 



(the semiclassical conductivity does not depend on q). 
Here r denotes the relaxation-time (RT), which in a 
phenomcnological way takes into account losses due to 
electron- impurity, electron-defect, and electron-phonon 
scattering. Equation ([5]) is obtained by assuming zero 
temperature T ~ 0, which is a good approximation for 
highly doped graphene considered here, since Ep 3> ksT. 
From Eqs. @ and ([5]) it is straightforward to obtain 
plasmon dispersion relation: 



q(u) = 
as well as losses, 



irh 2 e (e r i + e r2 ) 



(i + — w 

TLU 



$tq 
3<7 



27TCT 



(6) 



(7) 



In order to quantify losses one should estimate the re- 
laxation time r. If the frequency tu is below the in- 
terband threshold frequency Wi„t er , and if u < uiQph, 
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then both interband damping and plasmon decay via 
excitation of optical phonons together with an electron- 
hole pair are inactive. In this case, the relaxation time 
can be estimated from DC measurements [1, [li|, i.e., it 
can be identified with DC relaxation time which arises 
mainly from impurities (see Refs. UH). It is rea- 
sonable to expect that impurity related relaxation time 
will not display large frequency dependence. In order 
to gain insight into the losses by using this line of rea- 
soning let us assume that the doping level is given by 
Ep = 0.64 eV (corresponding to electron concentration 
of n = 3 x 10 13 cm -2 ); the relaxation time corresponds 
to DC mobility fj, = 10000 cm 2 /Vs measured in Ref. @: 
tdc = n%\/rvK I evp = 6.4 x 10~ 13 s. As an example, 
for the frequency tkx> = 0.155 cV (\ a ir = 8/im), the 
semiclassical model yields 'Slq/^sq w 151 for losses and 



42 for wave localization. Note that both of 



these numbers are quite favorable compared to conven- 
tional SPs [e.g., see Fig. [He)]. It will be shown in the 
sequel that for the doping value Ep = 0.64 eV this fre- 
quency is below the interband loss threshold, and it is evi- 
dently also smaller than the optical phonon loss threshold 
TvjJOph ~ 0.2 eV, so both of these loss mechanisms can 
indeed be neglected. 



Here f(E) = [ e {E-E F )/k B T + ^-l is the F errm distribu- 
tion function, Ep is the Fermi energy and factor 4 stands 
for 2 spin and 2 valley degeneracies. Note that in Eq. (|8|) 
uj is given an infinitesimally small imaginary part which 
leads to the famous Landau damping; that is, plasmons 
can decay by exciting an electron-hole pair (interband 
and intraband scattering) as illustrated in Fig. [21(c). The 
effects of other types of scattering (impurities, phonons) 
can be accounted for by using the relaxation-time r as 
a parameter within the RPA-RT approach [2j|, which 
takes into account conservation of local electron number. 
Within this approximation the 2D polarizability is 



(l + i/mT)x(q,uj + i/T) 
l + (i/uT)x(q,u + i/T)/x(q,0)' 



(10) 



The 2D dielectric function and conductivity are respec- 
tively given by (see [Uj]): 



e RP A{q,uj) = h — Xr{q,^), (H) 

z zeo 



and 



B. RPA and relaxation-time approximation 

In order to take the interband losses into account, we 
use the self-consistent linear response theory, also known 
as the random-phase approximation (RPA) [33[ , together 
with the relaxation-time (finite r) approximation intro- 
duced by Mermin (29[. Both of these approaches, that is, 
the collisionless RPA (r -> oo) [Mill], and the RPA-RT 
approximation (finite r) [20| , have been applied to study 
graphene. In the r — ► oo case, the RPA 2D polarizability 
of graphene is given by [22l |: 



(12) 



We note here that throughout the text only 7r-bands are 
taken into consideration; it is known that in graphite, 
higher er-bands give rise to a small background dielectric 
constant [33] at low energies, which is straightforward to 
implement in the formalism. Using Eqs. (0| and (TT2|) we 
obtain that the properties of plasmons (i.e., dispersion, 
wave localization and losses) can be calculated by solving 



ffiPi(?,w) = 0, 



(13) 



X(q,u)) = -Jl{q,u), (8) 



where 



with complex wave vector q = q\ +iq2- The calculation is 
simplified by linearizing Eq. (|13[) in terms of small 92/51, 
to obtain, 



. j_ ST^ /(^»2,k+q) - /(-5ni,k) 

Z-j fuj + E n k - £„ 2;k+q 
x |(m,k| e -^>2,k + q)| 2 . 



(9) 



e r i + e r 2 



2 2e QQl 
for the plasmon dispersion, and 



Rp(gi,w)]=0, 



(14) 



<Z2 = 



j3?[n( gi , W )] + ^[U( qi ,uj)(l - n(gi, W ))/n(gi,0)] 



ifkmiqi 



(15) 



yielding losses. Note that in the lowest order the disper- 
sion relation (and consequently X a ir/\ and the group 



velocity v g ) does not depend on r. This linearization 
is valid when q 2 <C q\, as the plasmon losses increase, 
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FIG. 3: (color online) Properties of plasmons in doped graphene. 
calculation, and the dashed lines with the semiclassical approach. 



Solid-lines are obtained with the number-conserving RPA 
Losses (a), field localization (wave "shrinkage") (b), and 



group velocity (c) for doping Ef = 0.135 eV, and relaxation time r = 1.35 x 10~ s, which corresponds to the mobility of 
10000 cm 2 /Vs. The upper scale in all figures is frequency v — uj/2iy, whereas the rose shaded areas denote the region of high 
interband losses. 



e.g., after entering the interband regime [the rose area 
in Fig. IHd)], results from Eqs. (HU) and (15]) should be 
regarded as only qualitative. The characteristic shape of 
the plasmon dispersion is shown in Fig. HJd). Note that 
the semi-classical model and the RPA model agree well if 
the system is sufficiently below the interband threshold 
[for small q, w(q) ~ y/q as in Eq. JBJ)]. By comparing 
Figs. [2jd) andQJb) we see that the dispersion for SPs 
on silver-dielectric surface qualitatively differs from the 
plasmon dispersion in graphene [34| . While SPs' disper- 
sion relation approaches an asymptote (to — > lusp) for 
large q values [Eq. Q], graphene plasmon relation gives 
u>(q) which continuously increases [Fig. HJd)]. 

Theoretically predicted plasmon losses ^Rq/^sq and 
wave localization A a i r /A p are illustrated in Fig. [3] for 
doping level Ep = 0.135 eV and relaxation time r = 
1.35 x 10 -13 s. We observe that for this particular dop- 
ing level, for wavelengths smaller than Xmter ~ 7.7 /zm, 
the system is in the regime of high interband losses (rose 
shaded region). Below the interband threshold, both 
losses and wave localization obtained by employing RPA- 
RT approach are quite well described by the previously 
obtained semiclassical formulae. Since the frequencies 
below the interband threshold are (for the assumed dop- 
ing level) also below the optical phonon frequency, the 
relaxation time can be estimated from DC measurements. 

At this point we also note that in all our calculations we 
have neglected the finite temperature effects, i.e., T w 0. 
To justify this, we note that for doping values utilized 
in this paper the Fermi energies are 0.135 eV~ 5.2kpT r 
(n = 1.35 x 10 12 cm" 2 ) and 0.64 eV« 25k B T r (« = 3x 
10 13 cm -2 ) for room temperature T r = 300 K. The effect 
of finite temperature is to slightly smear the sharpness 
of the interband threshold, but only in the vicinity (~ 
kpT r ) of the threshold. 

By increasing the doping, Ep increases, and the region 
of interband plasmonic losses moves towards higher fre- 
quencies (smaller wavelengths). However, by increasing 
the doping, the interband threshold frequency will even- 
tually become larger than graphenc's optical phonon fre- 



quency u>oph- there will exist an interval of frequencies, 
k>Oph < ^ < Winter, where it is kincmatically possible 
for the photon of frequency u> to excite an electron-hole 
pair together with emission of an optical phonon. This 
second order process can reduce the relaxation time esti- 
mated from DC measurements and should be taken into 
account, as we show in the following subsection. 



C. Losses due to optical phonons 

In what follows, we estimate and discuss the relaxation 
time due to the electron-phonon coupling. This can be 
done by using the Kubo formula which has been utilized 
in Ref. [3l[ to calculate the real part of the optical con- 
ductivity, ^Ra(uj,q = 0). The calculation of conductivity 
5icr(a;,0) involves the electron self-energy £(-E), whose 
imaginary part expresses the width of a state with en- 
ergy E, whereas the real part corresponds to the energy 
shift. Let us assume that the electron self-energy stems 
from the electron-phonon coupling and impurities, 

E(E) = Ze- P h(E) + E imp (E). (16) 

For T, e _ p h we utilize a simple yet fairly accurate model 
derived in Ref. [30|: If \E — Ep\ > hu> 0p h, then 

&E e - ph (E) = j\E - sgn(£ - Ep)faj ph\, (17) 

while elsewhere QY, e - v h(E) = 0; the dimensionless con- 
stant 7 = 18.3 x 10 -3 [301 i s proportional to the square 
of the electron-phonon matrix element f3Cj, i.e., the 
electron-phonon coupling coefficient. In order to mimic 
impurities, we will assume that ^~Ei mp (E) is a constant 
(whose value can be estimated from DC measurements). 
The real parts of the self-energies are calculated by em- 
ploying the Kramers-Kronig relations. In all our calcu- 
lations the cut-off energy is taken to be 8.4 eV, which 
corresponds to the cut-off wavevector k c — ir/a, where 
a = 2.46 A. By employing these self-energies we calculate 
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FIG. 4: (color online) (a) The real part of the conductivity in units of <ro = ire 2 /2h in dependence of frequency hu>/EF, and 
(b) the corresponding relaxation time as a function of wavelength. The contribution to R<t(oj) from impurities is chosen to 
be negligible. The displayed graphs correspond to two different values of doping which yield Ef — 0.135 eV (solid blue line), 
and Ef — 0.640 eV (dashed red line). The position of the optical phonon frequency huioph ~ 0.2 eV is depicted by the dotted 
vertical line in (b); dot-dashed lines depict the values of wavelengths corresponding to 2Ef, that is, the interband threshold 
value (for q — 0) for the two doping concentrations. 



the conductivity $ta(u>,q = 0), from which we estimate 
the relaxation time by using Eq. (0, i.e., 

I \ G^Ep 1 . . 

■nn z !jj z yta{tu, 0) 

for the region below the interband threshold; in deriving 
(UHJ) we have assumed rw> 1. 

Figure 2] plots the real part of the conductivity and the 
relaxation time for two values of doping: Ep = 0.135 eV 
(n = 1.35 x 10 12 cm" 2 , solid line) and E F = 0.64 eV 
(n = 3 x 10 13 cm -2 , dashed line). In order to isolate the 
influence of the electron-phonon coupling on the conduc- 
tivity and plasmon losses, the contribution from impuri- 
ties is assumed to be very small: ^sT, imp (E) = 10~ 6 eV. 
The real part of the conductivity has a universal value 
(To = ire 2 /2h above the interband threshold value htu = 
2E F (for q = 0), e.g., see [U HI & We dearly see 
that the relaxation time is not affected by the electron- 
phonon coupling for frequencies below uo v h, that is, we 
conclude that scattering from impurities and defects is 
a dominant decay mechanism for u < uJo P h (assuming 
we operate below the interband threshold). However, for 
to > uJoph, the relaxation times in Fig. 2]are on the order 
of 10~ 14 — 10~ 13 s, indicating that optical phonons are 
an important decay mechanism. 

It should be emphasized that the exact calculated val- 
ues should be taken with some reservation for the follow- 
ing reason: strictly speaking, one should calculate the re- 
laxation times t(u>, q) along the plasmon dispersion curve 
given by Eq. (| 14[> : namely the matrix elements which en- 
ter the calculation depend on <j, whereas the phase space 
available for the excitations also differ for q = and 
q > 0. Moreover, the exact value of the matrix clement 
for electron phonon coupling is still a matter of debate 
in the community (e.g., see Ref. [3f|). Therefore, the 
actual values for plasmon losses could be somewhat dif- 
ferent for to > woph- Nevertheless, fairly small values of 



relaxation times presented in Fig. @] for to > tooph indi- 
cate that emission of an optical phonon together with an 
electron-hole pair is an important decay mechanism in 
this regime. Precise calculations for q > and to > tooph 
are a topic for a future paper. 

Plasmonic losses and wave localization calculated from 
the RPA-RT approximation are illustrated in Fig. [5] for 
doping level Ep = 0.64 eV and the relaxation time r 
given by r _1 = t^q + T~} h , where t dc = 6.4 x 10~ 13 s 
(mobility 10000 cm 2 /Vs), whereas r e _p^ is frequency de- 
pendent and corresponds to electron-phonon coupling as- 
suming very clean samples [see dashed line in Fig. 0Jb)] . 
Interband losses [left (rose shaded) regions in all panels] 
are active for wavelengths smaller than Xmter ~ 1.7 //m. 
In the frequency interval u>%nter > w > UJoph [central (yel- 
low shaded) regions in all panels], the decay mechanism 
via electron phonon coupling determines the loss rate, 
i.e., t m T e -ph- For to < u>o P h [right (white) regions in 
all panels], the DC relaxation time tdc can be used to 
estimate plasmon losses. 

It should be noted that the mobility of 10000 cm 2 /Vs 
could be improved, likely even up to mobility 
100000 cm 2 /Vs 13], thereby further improving plasmon 
propagation lengths for frequencies below the optical 
phonon frequency. However, for these larger mobilities 
the calculation of losses should also include in more de- 
tails the frequency dependent contribution to the relax- 
ation time from acoustic phonons (this decay channel is 
open at all frequencies); such a calculation would not 
affect losses for to > u>oph where optical phonons are 
dominant. 



IV. CONCLUSION AND OUTLOOK 

In conclusion, we have used RPA and number- 
conserving relaxation-time approximation with experi- 
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FIG. 5: (color online) Properties of plasmons in doped graphene. Solid-lines are obtained with the number-conserving RPA 
calculation, and the dashed lines with the semiclassical approach. Losses (a), field localization (wave "shrinkage") (b), and 
group velocity (c) for doping Ef = 0.64 eV; losses are calculated by using the relaxation time r _1 = + T~_ ph , where 
TDC = 6.4 x 1CT 13 s, and r e _ p h is the relaxation time from the electron-phonon coupling for the given parameters. In the white 
regions (right regions in all panels), losses are determined by Tnc- ln the yellow shaded regions (central regions in all panels), 
losses are determined by the optical phonon emission, i.e., r e _ p h. The rose shaded areas (left region in all panels) denote the 
region of high interband losses. Dotted vertical lines correspond to the optical phonon frequency u>oph ~ 0.2 eV. The upper 
scale in all figures is frequency v — uj/2-k. See text for details. 



mentally available input parameters, and theoretical es- 
timates for the relaxation-time utilizing electron-phonon 
coupling, to study plasmons and their losses in doped 
graphene. We have shown that for sufficiently large dop- 
ing values high wave localization and low losses are simul- 
taneously possible for frequencies below that of the op- 
tical phonon branch ui < ujo P h (i-e., E p i asmon < 0.2 eV). 
For sufficiently large doping values, there is an interval of 
frequencies above uJoph and below interband threshold, 
where an important decay mechanism for plasmons is ex- 
citation of an electron-hole pair together with an optical 
phonon (for u> < LUoph this decay channel is inactive); 
the relaxation times for this channel were estimated and 
discussed. We point out that further more precise calcu- 
lations of plasmon relaxation times should include cou- 
pling to the substrate (e.g., coupling to surfacc-plasmon 
polaritons of the substrate), a more precise shape of the 
phonon dispersion curves [28[ , and dependence of the re- 
laxation time via electron-phonon coupling on q > (see 
subsection IHI C|) . 

The main results, shown in Figures [3] and O point 
out some intriguing opportunities offered by plasmons 
in graphene for the field of nano-photonics and meta- 
materials in infrared (i.e. for lo < u>oph)- For exam- 
ple, we can see in those figures that high field localiza- 
tion and enhancement X a ir/\ ~ 200 [see Figure [3Jb)] 
are possible (resulting in X p < 50 nm), while plasmons 
of this kind could have propagation loss-lengths as long 
as ~ 10Ap [see Fig. Ufa)]; these values (albeit at dif- 
ferent frequencies) are substantially more favorable than 
the corresponding values for conventional SPs, for exam- 
ple, for SPs at the Ag/Si interface X a ir/\ ~ 20, whereas 
propagation lengths are only ~ 0.1A sp [see Fig. [He)]. 



Another interesting feature of plasmons in graphene is 
that, similar to usual SP-systems fjj], wave localization is 
followed by a group velocity decrease; the group veloci- 
ties can be of the order v g = 10~ 3 — 10~ 2 c, and the group 
velocity can be low over a wide frequency range, as de- 
picted in Figs. OUc) and[^c). This is of interest for possi- 
ble implementation of novel nonlinear optical devices in 
graphene, since it is known that small group velocities 
can lead to savings in both the device length and the 
operational power [5(|; the latter would also be reduced 
because of the large transversal field localization of the 
plasmon modes. 
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